Gravitational collapse of charged dust cloud in the Lovelock gravity 
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We study the effect of charge on gravitational collapse of inhomogeneous dust cloud in the Einstein, 
Gauss-Bonnet and Lovelock gravity. Dynamics of the collapsing shell is analyzed. The conditions 
for the occurrence of bounce during collapse are given. We also show that shell crossing occurs 
inevitably if the shell near the center is weakly charged. 
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(N ! I- INTRODUCTION 

The final fate of gravitational collapse is one of the important unsolved problems in classical general relativity 
(GR). Both the possibilities, i.e., the formation of naked singularity or black hole as an end-state of collapse, may 
present a crucial step in our understanding of nature. On one hand, we have the singularity theorems which tell us 
that singularities are inevitably formed through gravitational collapse with physically reasonable and general initial 
data [![. These singularities lead to a loss of predictability if visible to a distant observer, and hence a break down of 
i— i- GR happens. On the other hand, the formation of naked singularities can imply an access to strong gravity regions. 
That is, they offer a possibility to probe the laws of nature and a window to a quantum gravity theory. However, one 
may wish that singularities should be wrapped by event horizons. This idea is paraphrased in the form of "cosmic 
censorship conjecture" (CCC) Q. If this conjecture is indeed true, we have black holes. A proof or otherwise of the 
, ^0 , CCC remains one of the unresolved issues in classical GR. 

The first pioneering work in gravitational collapse was carried out by Oppenhcimcr and Snyder Q . They considered 
7— I | the spherical collapse of a homogeneous dust cloud. They showed that black hole is formed and the singularity is 
^ ■ enclosed by the event horizon. However, many solutions which contain naked singularities have been found, for 
example, in the collapse of inhomogeneous dust cloud [il-fllj]. We now understand that, in the Oppenheimer-Snyder 
model, the homogeneity of dust cloud and spherical symmetry are important for the formation of black holes. In order 
to examine the generality of the above results, it is important to study the collapse with generic initial conditions. 
Then, gravitational collapse has been studied in various situations, e.g., collapse of charged dust cloud, dust with 
jn . tangential pressure, perfect fluid and so on [l2l - [r7j |. 

In the late stages of collapse with the strong gravity, GR does not provide us with an accurate description of the 
final fate of gravitational collapse. Motivated by this and the fact that the best candidate fundamental theory such as 
superstring/M theory lives naturally in more than four dimensions, considerable efforts have been made on studies of 
gravitational collapse in higher dimension al g ravitational theories. Gravitational collapse in the simple extensions of 
GR to higher dimension is studied as well jlSl — F20j | . It has been found that the singularity formed through the spherical 
collapse of inhomogeneous dust cloud cannot be naked in dimensions higher than five with smooth initial data. 
However, spacetime around the singularity is extremely curved and the reliability of the pure GR is doubtful. The 
physical motivation to go beyond GR comes from the fact that quantum gravitational effects / string corrections play 
a significant role in the region near singularity. Actually, superstring theory predicts certain higher order curvature 
corrections to GR. From this fact it is important to study the gravitational collapse in theories with higher curvature 
corrections. 

In this paper we focus on a special type of gravitational theory with higher curvature corrections, that is the Lovelock 
theory [2l| . The Lovelock theory is the most general gravitational theory in which: (a) Lagrangian is covariant and 
(b) field equations contain up to second order derivatives of metric tensor. The Lovelock gravity is a widely studied 
higher dimensional gravitational theory. The black hole solutions were found in Refs [22H25l |. Some authors have 
studied the dust collapse in a special case of the Lovelock theory (2(| [53] ■ In our previous paper [28[ , we studied 
the spherical collapse of dust cloud in the full Lovelock theory with any spacetime dimension. We showed that the 
singularities are always formed and can be naked for some initial data. In addition we showed that the nature of 
singularity depends on the dimensionality of spacetime, that is, whether it is even or odd. 

As mentioned before, it is important to investigate the collapse in more general situations from the CCC point 
of view. To this end we will consider the spherical collapse of a dust cloud coupled with the Maxwell field, i.e., a 
charged dust cloud, in the Einstein, Gauss-Bonnet and Lovelock gravity. The collapse of a charged dust cloud can 
be regarded as a toy model for understanding rotating spacetimes. We will also compare the dynamics of a charged 
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dust cloud with a neutral dust cloud. In the charged dust case, using a kinematical argument, we see that the shell 
near the center inevitably bounces if the charge is sufficiently small. We also show that the shell-crossing singularities 
occurs before the shell bounces. This feature is quite different from the neutral dust cloud cases 28]. Although we 
have included higher order curvature corrections, our result is same as in the case of the Einstein gravity(See Refs. 
[29l |30|.). This is quite impressive because the higher curvature corrections do affect the dynamics in the case of 
neutral dust cloud. 

The rest of this paper is organized as follows. In Sec. II, we briefly review the Lovelock theory. In Sec. Ill, we derive 
the basic equations and solve these partially to discuss the shell motion kinematically. Looking at the key equations, 
in Sec. IV, we show that the shell-bounce will occur kinematically. We investigate the Einstein, Gauss-Bonnet and 
full Lovelock gravity case in turn. In Sec. V, we show that shell-crossings inevitably occur before the occurrence of the 
shell-bounce. In Sec. VI, we briefly summarize and discuss our result. In the appendices, we consider the matching 
condition of the inner collapsing with the outer electro- vacuum region and a detail of the calculation in the text. We 
use the geometrized units, G = c = 1. 



II. LOVELOCK GRAVITY: THE EXTERIOR GEOMETRY 



In this section, we briefly review the Lovelock theory of gravity in D = (n + 2)-dimensional spacetimes. The 
Lagrangian of the theory is 
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C=J2— An, (1) 
m—l 

where 

Here R pq sr is the Riemann curvature tensor of spacetime and S^i^'.'-i^^-i^Tm * s ^ ne generalized totally antisymmetric 
Kronecker delta. a m are arbitrary constants which cannot be determined by the theory. The suffixes /ii , . . . , ^ 2m and 
Ui,..., V2m run from 1 to D, and k is given by k = [(D — l)/2] ([x] is the integer part of x). 
The field equations derived from this Lagrangian are of the form 

fe 

C, v — _ \ * m XW\V2 — V3m-\V2m E> MlM2 D M2rn-lM2m 

A* 2 m+1 m WlM2 --- M2m - lM2m 2 ' ' ' V2m ~ ll/2m 

m—l 

= s^t;. (3) 

Here is the energy-momentum tensor of matter fields. 

The static black hole solution with spherically symmetry was found in Refs. [23Tj25j for vacuum and in Refs. (3ll.l32j 
for electro-vacuum cases. An electro-vacuum spacetime is described by the metric of the form 

ds 2 = ~f(r)dt 2 + -r^-dr 2 + r 2 f ii dx i dx j . (4) 
f(r) 

Here 7ij is the line element of n-dimensional maximally symmetric surfaces and / is 

f(r) = K~r 2 i>{r). (5) 
The curvature constant k can take values —1,0 and 1. The function ip(r) is a solution of the algebraic equation 

■ 2m-2 



m=2 L p=l - 1 



where /1 and q are constants proportional to the Arnowitt-Deser-Misner (ADM) mass and charge, respectively. This 
solution describes the spacetime exterior to a charged dust cloud (See Appendix A for the details). In what follows 
we shall restrict ourselves to the case k = 1 because we are interested in asymptotically flat spacetimes here. The 
generalization of our study to other cases with k = 0, — 1 is rather easy. 
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III. BASIC EQUATIONS: THE INTERIOR GEOMETRY 

In this section we shall derive equations for the interior of a collapsing charged dust cloud. The spacetime metric 
under the spherical symmetry ansatz can be written as 

ds 2 = -A 2 dt 2 + B 2 dr 2 + R 2 dQ 2 n , (7) 

where the metric functions A(t,r),B(t,r) and R(t,r) are arbitrary functions of coordinates t and r, and df^ is the 
line element of a n-sphere. Using the gauge freedom we set i?(0, r) = r on an initial data surface. We also assume 
that the interior of the collapsing cloud is composed of charged dust, i.e. the energy momentum tensor T^ v takes a 
form 

Here T^ D ^ V is the energy momentum tensor for a dust cloud 

T< D K = e{t,r)ui i u v , (9) 

and T^» v is that of the Maxwell field, 



(10) 



e(t, r) is the energy density of the dust cloud, is its 4- velocity and F^ v is the field strength of the Maxwell field. 
In addition, we have the Maxwell's equations given by 

Vj.F"" - -A-Kf , (11) 

and 

VfrKp] = • (12) 

Here f is the electric current and is the covariant derivative with respect to the spacetime metric. We assume 
that the current has the form 

r = ™ v , (13) 

where a is the electric charge density. 

In the case of a non-zero current, the Maxwell's equations become 



ABR 



1 Att 

' d r (ABR n F tr ) = —a (14) 



and 



We can solve these equations as 



ABR 



! d t (ABR n F tr ) = 0. (15) 



n Q(r) 



2 ABR n ' 

where Q(r) is the integration function and it is the charge contained in the ball of radius r, and it satisfies 



(16) 



Q' = —Ana BR". (17) 
n 



The prime stands for the differentiation with respect to r. 

The conservation equations of the energy momentum tensor, V„T^ = 0, lead to 



V M (eu") = (18) 



and 

eti"V„u /1 = oF» v u v . 
Using Eqs. jJJ), (|16j) and (fTTj) . the conservation equations become 

d t (eBR n ) = 

and 

'A'\ (n\ 2 QQ> 



e 



Eq. (|20|) can be integrated as 



A y V 2 / 4ttR 2 



N'(r) = -AirsBIF, 
n 



where N{r) is an arbitrary function of r. Using Eq. (|22p . we can re-express Eq. (j2"Tj) as 

A' _ nQQ N B 
~A ~ 2 i?" ' 

where 

Q!_ = °_ 
N' ~ e' 

Now, let us consider the field equations (|3]). They have the following nontrivial components: 

/ n \ TTl—X 

ST ^ ( R\ ( R '\ 2 1 

m—1 p—0 



Qn 



AR \BR R 2 



2 



X 



R_ AR A'R' (n- (2m- 1)) / / R V f R' \ 2 l_ 
A 2 R ~ A 3 R ~ AB 2 R + 2m [ AR j ~ { BR J + R 2 I J ~ \ 2 J R 2r 



f2 \ m—1 



m—1 p—0 
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R" B'R' BR (n - (2m - 1)) / / R \ ( R! \ 2 J_ 
B 2 R + B 3 R + A 2 BR + 2m 1 AR } \BRj + R 2 



87re(i, r) 



and 



k 2m-2 




m—1 p—0 

In the above, the dot stands for the differentiation with respect to t. 
First we realize that, for a non-trivial solution, Eq. (|27p implies 



73 2 i? AB 2 R B 3 R 
Using Eq. ([23j), we can easily integrate Eq. (l28t as 
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where W{r) is the integration function which corresponds to the total energy, as we will show later. 
Now we introduce a useful function L defined as 
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2m 11 v ^ U V 2(n-l)R n ~ 1 J R 2r 

m=l p=0 \ \ / - / 

Using this definition of L, we can rewrite Eqs. (|25|) and (|26|) in a simple form, i.e. 

1 d , Tm+ n_fn\ 2 Q 2 



RR™ dt x ' V 2 / R 2n 

and 

i a^,„-M, . /«r ^ 

R'R n dr 

The integration of Eq. (f3"Tj) gives us 



^ n+i ) = U ^77 (31) 



^" +1 )=8-+Q) ^. (32) 



V2J (n - l)i? 2 " 2 JfH-i ' [66) 

where F(r) is the integration function which corresponds to the mass, as shown later. Substituting this into Eq. (|32l) . 
we see that 

n Fjr)' n 2 QQ' 

2 2(n- l)i? 2 «-!i?'' 1 ; 

From this expression, we see that the spacetime admits two type of singularities in general. The shell-focusing 
singularity at R = and the shell-crossing singularity at R' — 0. Using Eq. (|29[) . we can re-express mass the function 
F(r) as 

F(r)' = 2WN(r)'. (35) 
Combining Eq. (|30|) with Eq. ([33| . we can obtain the basic equation which determines i? 

' 1 (,„ nQQ N \ 2 \ 1 /n\ Q 2 F{r) 



771—1 

where q is defined by 

- J 1 if 1 = 1 m 

l ~\fUl^in-p) if 2 < Z < A:. (37) 

In order for the initial surface to be regular at the center (r — 0), the functions F(r), Q{r) and W(r) should behave 
as 

F{r) =r n+1 M(r), (38) 
Q(r) = r n+1 q{r) (39) 

and 

W{r) = 1 + r 2 b(r), (40) 

where M(r),q(r) and b(r) are regular functions on the initial hypersurface, and it is supposed that q(r) and M(r) 
behave as q'(r),M'(r) = 0(r) near r = 0. 

Before closing this section, we would like to make a remark. One might wonder whether 

k 2m-2 / / p \ / R' \ 2 1 \ 

S-n(-rt((^)-(^)^] -0 
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is also a solution to Eq. ([27]) . But, combining with Eq. (|25|) . we have 

II («-*) ^ 

m— 1 p— 

We can see that Eqs. (|26[) and (l4"2"j) imply e(t, r) — 0. Hence, the spacetime must be electro-vacuum which is of no 
interest to us here. 




IV. ANALYSIS OF SHELL MOTION 



In this section we shall study a weakly charged shell in the vicinity of the center (r = 0) and show that its bounce 
is inevitable. The analysis in this section is also used to show the occurrence of shell-crossing singularity in the next 
section. For a pedagogical reason, we shall start with the Einstein gravity case and then discuss the Gauss-Bonnet 
and full Lovelock gravity cases. 



A. Einstein gravity case 



Here we investigate the spherical collapse of charged dust in the Einstein gravity with arbitrary dimensions. In this 
case, Eq. ([55)1 becomes 



= -V E , 



(43) 



where 



V E = -E(r) - 



M(r) 
R"- 1 



Q 2 



(n- l) 2 



Qn 



2(n- 1) 



R 2(n-1 



In the above E{r) = W 2 {r) - 1 and M(r) = F(r) - nWQQ N /(n - 1). 

Let us consider whether shell bounce (Ve = 0) occurs or not. From the discriminant of Ve 
we find that there is no solution of R for Ve = if 



(44) 



with respect to R, 



M 2 {r) < 



4E(r)Q 2 (r) 
(n - l) 2 



Qn 



2(n- 1) 



(45) 



Thus, if Eq. (1431) is satisfied, (i) shell always collapses into shell-focusing singularity(i? = 0) because Ve stays negative 
definite throughout collapse in the E > case, and (ii) there is no solution for Eq. (|4"3")l in the E < cases because 
Ve always has a positive value. 

Similarly we can analyze different cases for the condition 



M 2 {r) > 



AE(r)Q 2 (r) 
(n- l) 2 



Qn 



2(n- 1) 



(46) 



In this case there are various possible situations depending on the sign of E, M and Q 2 N — 2(1 — 1/n). For example, 
in the case for E > 0, M > and Q 2 N — 2(1 — 1/") > 0, Ve always has negative value. This implies that the shell 
can hit the shell focusing singularity. On the other hand, in the case for E > 0, M > and Q 2 N — 2(1 — 1/n) < 0, Ve 
cannot have negative value for small R. This means shell inevitably bounces. In Table 1, we summarize the result 
satisfying the condition above. 





E > 


E = 


E < 




M > 


M < 


M > 


M < 


M > 


M < 


n i ^ 2(n-i) 

VjV > n 


singular 


singular and bounce 


singular 


singular 


singular 


singular 


n 2 „ 2(n-l) 


bounce 


bounce 


bounce 


no solution 


oscillations 


no solution 



Table 1 The cases satisfying condition Eq. (|46|) . 
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In summary we could show that (i) every shell with Q N < 2(1 — 1/n) bounces if the solution to Eq. (|43|) exists, and 

2 



(ii) every shell with Q N > 2(1 — 1/n) collapses into shell- focusing singularity. In the Einstein gravity, the condition, 
Q 2 N < 2(1 — 1/n) seems to play a significant role for occurrence of the bounce, 
condition as 



From Eq. (f22[) . we can express this 



1 

1 

n 



1/2 



(47) 



This means that a weakly charged collapsing shell will inevitably bounce in the case of the Einstein gravity. The 
condition for the bounce depends on the ratio of the charge to the dust density and the number of spacetime dimensions. 
Note that this result is basically same as that obtained for four dimensions [29|, |3(| • 




B. Gauss-Bonnet/Lovelock gravity 

Next we shall consider the full Lovelock gravity cases including the Gauss-Bonnet one. Eq. becomes 

)m 
R2m \ 2 ) ( n -l)R 2n V ' 

Let us consider a certain epoch t = when 

i?™- 1 = V*L (49} 

2{n-l)F [ ' 

holds. At this epoch, the right-hand side of Eq. (|48|) vanishes. To achieve this condition in collapsing situations, the 
condition R n ~ l > 2 (n-i)F s h° u ld be satisfied at the initial surface. We can see that this is always satisfied near the 
center (r = 0) because the initial condition for R is R = r and the regularity condition for F(r) and Q(r) implies that 

2("-i)F S oes as C( r " +1 ) near tne center. 
At t = we can rewrite Eq. (|4*51) as 

FQ N \ 2 \ I 




rn—1 

Here we focus on the weakly charged cases which satisfy the condition 



2W F ~ Q ~ 2W F 

and one can show 



W-IF'Q'W + IF' 

< 77 < -TTF77 — 77 - ( 51 ) 




(52) 

But this is impossible. Therefore we may conclude that the shells near the center will inevitably bounce at least 
before t = t*, if the artificial condition (ISTl) is satisfied. 

Next, let us briefly examine the physical interpretation of the condition ([51]) . In the vicinity of the center, W(r) 
behaves as W(r) = 1 + b(r)r 2 and then the condition (fSTj) can be written approximately as 

O' F' 

0<|< T . (53) 

Roughly speaking, this is satisfied when the spatial scale of the mass distribution is shorter than charge one. The 
integration of the above gives us 



Q(r) < g(0) = 1 a 
F(r) ~ F(0) 2 e 



(54) 

We see that, in order to show the occurrence of the shell bounce, we needed somewhat artificial condition of Eq. 
(1511) which is not required in the Einstein case. We may conclude that the shell near the center bounces inevitably 
for sufficiently small charge in the sense of the condition (1541 . 
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SHELL-CROSSING 



In the previous section, using the kinematical argument, we showed that the shell bounces if the conditions (j47J) 
and (fSTj) are satisfied for the Einstein, Gauss-Bonnet and Lovelock gravity respectively. Hereafter, we assume these 
conditions. In this section we will show that the shell-crossing occurs near the center before the shell bounces if the 
charge is relatively small. In order to show this, we shall adopt a new coordinate system, called Mass- Area coordinates. 
This helps us to study on the shell crossings as in the four dimensional case of charged dust collapse [2^, [3(| • 



A. Mass-area coordinate 

We first employ the coordinate transformation (t, r) — > (F, R), where F is the mass function introduced in Eq. ([33 
and R is the area radius. Then, under the above coordinate transformation, the metric transforms as 



ds 2 = -A 2 dt 2 + B 2 dr 2 + R 2 dVL r , 
B ,2 



B 



dR 2 



B 



-A*f R f F +[ — ) [l-Rf R )Rf F 



dRdF + R 2 dfl 2 „ 



dF 2 



g RR dR z + 2g RF dRdF + g FF dF z + R 2 dn 2 



(55) 



where t = f(R, F), r = r(R, F), f R = d R f and f F = d F f. The reversibility of the coordinate transformations gives us 



JrR - l, 

r R = 0, 
F'r F = 1 



and 



f F F' + f R R' = 0, 

where r R = d R r and r F — d F r. The metric components can be written as 

g FF = -T 2 (u 2 (u R ) 2 ) , 



9RF 



and 



1 



9RR 



(u R ) 2 ' 

where u R = R/A is the i?-component of the shell velocity, and u and T are defined by 

nQQ N 



Rf 

u = — = W - 

B 2{n-\)R n - 1 



(56) 
(57) 
(58) 



(59) 

(60) 
(61) 

(62) 

(63) 
(64) 



Using Eq. (fT9)) . we can show 



Fr, 



u{u R ) 2 ' 



(65) 



where F R = d R F. See Appendix B for the details of the derivation. 
Let us define ip by 



{u R ) 2 + 1 - u 2 
R 2 



(66) 
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Note that V> satisfies the following equation which is same as Eq. ([36 



m= 1 



Differentiating Eq. ([55)1 with respect to F, u R can be written as 



u, 



iti? is computed from the definition of u (Eq. (|63l) ) as 

^ = ^- 4(n-l)M/fl"-i ' (69) 
where we have used dN/dF — 1/2W, which is directly derived from Eq. (j3"5")l . Now, we can rewrite Fr as 

r* = -T^(wr- Sm^ + ^)- <™> 



R) 3 V 4(n- 2u 



(« J 

Then the integration of Eq. (|70|) gives us 

where % is the integration function of F. 

From Eq. (|59l) we can see that T = corresponds to i?' = i.e. shell-crossing singularity. In the remaining part, 
we investigate whether the shell-crossings occurs or not for each cases. 

Note that the energy density in the current coordinates is written as 

n 1 

£ ~ ~ 4 AitW R n u R J- ^ ' 

We observe that T should be positive in collapsing situations because it is assumed that u R is negative and e is 
positive. From now on we consider the marginally bound case (W = 1) for simplicity. 



B. Einstein case 



In the Einstein gravity case, ip is 



and ipF becomes 



F nQ 2 
ii^+i ~ 2(n- 1)R 2 



V- = T^TT ~ — TTTPIn ( 73 ) 



R n+1 (n - l)R 2n R^ 1 
where we used the fact that Qf — Qn/2 holds for W = l(See Eq. ([35]) ) and the definition of u(Eq. ([63]) ). 



Then Eq. (fTTj) can be re-written as 

dR 



2{u R fR n - 1 



(l - 2{n n _ 1} {Q% + QQnn)^ + H{F). (75) 



Here we note that = O(r ) and QQwat = C(r 2 ) because of Q = 0(r n+1 ) and N = 0{r n+1 ) (see Eqs. (35J), (05]) 
and ([39])). Thus, if the condition 



(J 



2 



< 2(1 - l/n) (76) 



is satisfied, the integrand is positive near the center. 

Since u R = ^ ~ near the bounce, T diverges with the negative sign. We see from Eq. (|72l that T is positive 
at the initial surface, but _F is negative near the shell bounce from the above argument. This implies that the shell 
crossing singularity occurs near the center before the shell bounces if the charged cloud is weakly charged, in the sense 
of |f | < (2-2/n) 1 / 2 (Eq. (SID). 
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Gauss-Bonnet case 



In the Gauss-Bonnet case, ip is 



12 \ 1/2 



1 1 / AaF 2anQ 2 , 
^ = -2^ ± 2^{ 1+ R^- (n-l)R^ ) ■ (77) 

From the condition (f5T|) . ip is restricted to positive values for the shell near the center l . Thus we only consider the 
plus branch of Eq. (J77J). 

From the direct computation, we can see that 

, 1 (, n QQn\ 

^ F = '/ 7TT72 ( 1 - ofa _ ia on-i ) ( 78 ) 

K \ L + BFFT - („-l)JR: 

and then Eq. (jTTj) becomes 

F -l 21^ 1 ; ~ ' ^ -^T) C«r + <W M ) 1 + H(n (79) 




(n-l)_R 2 



We observe that if the condition 



— ~ 9 < — 1 wo (80) 

£ « fit 4aF _ 2anQ2 \ ' H - 1 A 2a(n-l)FN 1/2 



is satisfied, the integrand of Eq. (|79|) will be positive near the center. Here, R* is defined as i?™ -1 = („_i"(„ + i) 

which maximizes the term 1 + jj0k ~ (n-i)R in ^ or gi ven F an d Q- This leads to the divergence of the integral in Eq. 
([?9"f with the negative sign. We saw from Eq. ([721) that T is positive at the initial surface. On the other hand, T will 
be negative near the shell bounce. Therefore, we can conclude that the shell near the center hits the shell-crossing 
singularity before the shell bounces if the conditions (f5T|) and ([50")) are satisfied. 

Now let us examine the meaning of the conditions (ppTj) and ([50)1 again. First, the condition (150"f corresponds to 
Q 2 N < 2 ^ n ~ 1>> in the Einstein gravity. This means that the shell is weakly charged compared to mass (see Eq. (|47j) L 
Next, under the condition (|80[) . the condition (|5Tj) roughly implies that the charge contained in the ball of radius r is 
smaller than the mass. This can be seen from Eq. ([54)) . As explained in Sec. IV, we need this condition to show the 
occurrence of shell crossing singularity. 

In summary, we can conclude that the shell near the center (r = 0) hits the shell-crossing singularity if the conditions 
(|5Tj) and ((80j) are satisfied. 



D. Lovelock gravity case 



In the Lovelock gravity case, T is given by 

R dR 



and where 



1 = - I U^rnLr- 1 " 2(^1) ^ + QQNN ^) + (81) 



, = 1 L _ n QQ N . 

VF R n + 1 ^ =1 mc m i> m - 1 V 2(n - 1) i?"- 1 ' ' 



1 First we can see that QQm — r n+1 c 2 (0)/e(0) is positive near the center (r = 0). And it is easy to see that the function Z := 
1 — (1 — 2 (n— 'l^"- 1 ) 2 ^ as on ^y one max i mum value with respect to _R, and it tends to zero as R — > oo and diverges with a negative sign 

as R — > +0. In addition, at t = t*, Z has a positive value of 1 — (1 — 2Q' F/QF') 2 because of the condition l|51|l . Since ip = jgzi^s + Z), 
it turns out that ip is positive from the beginning to the epoch t = t„ near the center. 
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This is obtained by differentiating Eq. (|47l) with respect to F. Because of the condition (|51l) . ip is restricted to positive 
values for the shell near the center (See the footnote around Eq. (1771) ). As in the Gauss-Bonnet case, if the condition 



is satisfied, the integrand of Eq. (|8Tj) will be positive near the center. From Eq. (|72|) . J- is positive at the initial 
surface, but we saw that F will be negative near the shell bounce. Therefore we can conclude that the shell near the 
center hits the shell-crossing singularity if the conditions (f5Tj) and are satisfied. 

The meaning of the conditions (1ST1) and (|83l) is the same as in the Gauss-Bonnet case. The condition (|83[) corresponds 
to Q 2 N < 2(n—\)/n in the Einstein case, which means dust is weakly charged as compared to mass. Then the conditions 
(|5T|) and (|83|) roughly imply that the charge contained in the ball of radius r is smaller than the mass. 



In this paper we considered the spherical collapse of a charged dust cloud in the Einstein, Gauss-Bonnet and 
Lovelock gravity. In the case of spherical collapse of a dust cloud [28[ , all shells of the dust cloud cannot bounce and 
will hit the singularity. But, in the case of collapse of a charged dust cloud, we found that the weakly charged shells 
near the center can bounce kinematically in each theories of gravity. This is because of the existence of a repulsive 
force exerted by the charge. We also found that such shells inevitably hit the shell focusing singularity in the theories 
of gravity considered here. 

In the higher dimensional Einstein gravity case, we found that the the shells near the center satisfying the condition 
Q 2 N = (|-) < 2(ro — l)/n, where a is the electric charge density and e is the mass density, can bounce kinematically. 
We also found that such shells inevitably hit the shell crossing singularity. In the four dimensional case, it has been 
shown that the shells with sufficiently small charge inevitably bounce [ID, [3(| • We found that this feature remains to 
be correct even in the higher dimensional Einstein gravity. 

Next we considered the Gauss-Bonnet and Lovelock gravity cases. In both of the theories, we found that the shell 
near the center can bounce under the condition ()51[) . which is not imposed in the Einstein gravity. We also found 
that the shell hits the shell-crossing singularity with the condition (j5"Tj) and ([50]) in the Gauss-Bonnet case, (|5Tj) and 
(|83f in the Lovelock gravity case respectively. We examined the meaning of the conditions (ISTl) . (|80j) and (f83|) . and 
it turns out that these conditions can be interpreted as that the shell is weakly charged compared to mass. Then we 
can conclude that the feature of the shell dynamics in four dimension is shared even in the Lovelock gravity case. 

We note that the reason why we imposed the additional condition (ISTj) in the Gauss-Bonnet and Lovelock gravity 
cases is purely technical. Therefore, our results should be true for more general initial data. The condition should 
be relaxed if we choose the epoch more sophisticated than t = t* which we chose in this paper, or if we perform a 
numerical analysis. The above remaining issue will be addressed in the near future. Since our initial data are specific, 
it would be interesting to see if there exists initial data sets which avoid shell crossing singularities in the Gauss 
Bonnet and Lovelock gravity theories analogous to the four dimensional results [13] • 
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In this appendix, we will address the junction between the inner and outer solutions. We assume that the inner 
solution is the collapsin g du st cloud presented here and the outer solution is the static vacuum black hole. We follow 
the argument in Refs. [26l |33| where the junction conditions in the Gauss-Bonnet theory and Einstein theory are 
discussed. 

Let £ to be the boundary of the two regions. The metric of the inner solution is 




(83) 



VI. SUMMARY AND DISCUSSION 
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Appendix A: Junction between inner and outer solutions 



ds 2 = ~A 2 dt 2 + B 2 dr 2 + R 2 da 



'71 



(Al) 



12 



and that of the outside region is 



ds 2 = -f{R)dT 2 + -j^rdR 2 + R 2 dfl n . (A2) 

i{R) 



We suppose that £ is described by the parametric equations R — R^(t) and T = In addition, it is natural to 

think that the boundary is comoving, that is, r = =constant. The induced metric q on £ from the metric of the 
inner region is written as 

q = -A 2 dt 2 + R^dfl n . (A3) 
On the other hand, using the metric of the outer region, it can also be rewritten as 



2 



q = -\f(R^)f 2 -^-\dt 2 + Rldn n . (A4) 



Since these should be identical, we see that 



A ° = [ nR ^-mr)l (A5) 



holds. 

The extrinsic curvatures of £ evaluated from the inner metric are 



AB 
and 

R' 



Kl = -4s (A6) 



-K} = ^. (AT) 
In terms of the metric of the outer region, the extrinsic curvatures also have another expression 

1 / af . . f 3 \ 

+K* = : ^ [RT-TR+ -±jR 2 T + —ff (A8) 



HI 



(f(Rz)T 2 f(Rs) 

and 



3 /2 I 2f 2 



The continuity of the extrinsic curvature implies 



- (A9) 



and 



•/ A ' ■ B . 
R' = —R+—R' 
A B 



■ \ 2 



R\ (R! 



A 



B 



(A10) 



(All) 



where 



Then, from Eq. (IA11[) . we see 



/(ife) = l-iZf^(ife). 



(A12) 



R 



HA) + ^-(^) = ^ 



(A13) 
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We recall that -0 satisfies Eq. (j67j 



and R is the solution to Eq. ([36 




2 \ 2m 



^ , , i? \ 1 (R'\\ _F(r ) Q(r ) 2 



It can be easily checked that Eq. (|Al3j) holds only if 

/x = F(r ) (A16) 
q = Q(r ). (A17) 

This shows that the inner region can be naturally joined to the outer region if ^ = F(tq) and q — Q(tq) is satisfied. 

Appendix B: Equation of motion of charged shell in the mass-area coordinate 

In this section, writing down Eq. p^|) in the mass-area coordinate (M, R) , we derive Eq. {SS"]) . We first note that 
u R = R/A and the other components of vanish. The non-trivial component of Eq. (fH?)) is 



Now, since 



and 



eu^Vf.u" = aF R R u R . (Bl) 



F\ = Rr R F\ + R'f R F^=l-^ (B2) 



/,' uR R 



^RR 



" (7^ + ^^^^). (B3) 



l r Fu R \{ 



eu ( u R F ,\ nuQ 



Eq. (IBip gives us 

y((# + ( ^ + ^ = CT iS- (B4) 

Using u, R = nQQ N /2R n (See Eq. O) and Q N = cr/e(See Eq. {23])), we obtain Eq. {55]) 
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